The Massive Schwinger Model 
a Hamiltonian Lattice Study in a 
Fast Moving Frame 

Helmut Kroger and Norbert Scheu 

Departement de Physique, Universite Laval, 
Quebec, Quebec G1K 7P4, Canada, 
Email: hkroger@phy.ulaval.ca, nscheu@phy.ulaval.ca 

February 7, 2008 

Abstract 

We present a non-perturbative study of the massive Schwinger 
model. We use a Hamiltonian approach, based on a momentum lat- 
tice corresponding to a fast moving reference frame, and equal time 
quantization. We present numerical results for the mass spectrum of 
the vector and scalar particle. We find good agreement with chiral 
perturbation theory in the strong coupling regime and also with other 
non-perturbative studies (Hamer et al., Mo and Perry) in the non- 
relativistic regime. The most important new result is the study of the 
^-action, and computation of vector and scalar masses as a function 
of the 0-angle. We find excellent agreement with chiral perturbation 
theory. Finally, we give results for the distribution functions. We com- 
pare our results with Bergknoff's variational study from the infinite 
momentum frame in the chiral region. 







1 Introduction 



In order to study QCD non-perturbatively, lattice gauge theory has been 
most successful. However, there are some observables, where computational 
progress in lattice gauge theory has been slow. Examples are: Higher excited 
states of hadrons and mesons, finite density thermodynamics (quark-gloun 
plasma), dynamical scattering calculations of cross sections or phase-shifts 
and hadron structure functions, in particular in the region of small Q 2 and 
small xb (10~ 2 to 10 -5 ). In our opinion a non-perturbative Hamiltonian 
approach may be a viable alternative. In deep inelastic lepton-hadron scat- 
tering, the success of the parton model suggests the physical idea to use a fast 
moving frame also for a computational study of those processes. The parton 
model can be justified using the operator product expansion. In equal-time 
quantization, the Breit-frame is the most convenient choice of frame in order 
to interpret structure functions. In Ref.JI| the authors have proposed such 
a scheme based on equal-time quantization, using a lattice Hamiltonian on 
a momentum lattice corresponding to a fast moving frame (Breit-frame). It 
has been applied to the scalar 4 theory in 3 + 1 dimensions. Distribution 
functions and the mass spectrum in the close neighbourhood of the critical 
line of the second order phase transition have been computed. 

Here we study QED\ + i, the so-called massive Schwinger model, in this frame- 
work. This model is physically interesting due to the following properties: 
(i) It is super-renormalizable. (ii) It has a limit where the model is analyti- 
cally soluble: the massless limit [[J. (iii) It has an (axial) anomaly, (iv) One 
can study the topological properties associated with ^-actions and #-vacua. 
Most of the analytic work on this model has been presented in Refs.||-[[J. 
This model is a widely studied benchmark problem for new methods. It has 
been studied by a variational method in the infinite-momentum frame [[K] 



in discretized light-cone quantization [12], on the light-cone using orthogonal 



functions [IB], on a finite lattice in equal time quantization [TL\ and by strong 



coupling series expansion |14 



The purpose of this work is to show: (a) The use of a fast moving frame, 
such that v < c, in conjunction with equal-time quantization works well also 
in QEDi + i. It is not necessary to go to the infinite momentum frame or 
to quantize on the light-cone. We obtain quite precise results in the ultra- 
relativistic region m/g — > 0. It has been claimed that low energy states of 
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the massive Schwinger model are not well aproximated in equal-time field 
theory. We show here that the problem is not equal-time quantization, but 
the use of the rest-frame, which we avoid to do here, (b) We consider as 
most important new results of this work the non-perturbative results of the 
dependence of the vector and scalar mass on the #-angle. 

2 Method 

Starting from the Lagrangian, we use the axial gauge A 3 = to obtain the 
Hamiltonian 

H = [ L dx 3 {^ 3 id^ + + C [ L ^ 3 (^)-^f(V>V)- (1) 

J-L 2 J~L —0*3 

One introduces a space-time lattice given by spacing a with N = — lattice 
nodes. Via discrete Fourier transformation one goes over to a momentum 
lattice, with cut-off A = - and resolution Ap = j. Motivated by the parton 
picture we make the assumption that left-moving particles are not dynam- 
ically important, if physical particles are considered from a reference frame 
characterized by a velocity v — which is close to the velocity of light. 
We thus consider a momentum lattice where < p°,p 3 < P. In order to 
minimize the number of virtual particle pairs created from the vacuum, we 
choose a small lattice size. The reason for this is that the number of vacuum 
pairs is roughly proportional to a vacuum density times the lattice size. On 
the other hand, a fast-moving physical particle is Lorentz contracted; thus it 
fits in a small lattice volume (when compared to the rest-frame). 

For the purpose of computing the mass spectrum, we need to determine 
the vacuum energy. Because the vacuum has the quantum number P = 0, 
the vacuum energy (and only this) is computed in the rest frame. Because the 
model is super-renormalizable, one can perform the continuum limit a — > 0. 
The only renormalization necessary is the subtraction of the vacuum energy. 
On a space-time lattice, one has to satisfy a physical condition (scaling win- 
dow) a < £a < L, where £ is the correlation length in dimensionless units, 
related to the physical mass of the ground state by M = ^. In a strongly 
relativistic system, McP, thus when a — > the scaling window is replaced 
by XjP < 1/M < L. For more details compare with Ref. [0. 
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3 Numerical results 



3.1 Mass spectrum 

We diagonalize the Hamiltonian in a sector with momentum P = to obtain 
the vacuum energy E vac . Then we diagonalize the Hamiltonian in a sector 
P^O corresponding to a relativistic velocity. This yields an energy spectrum 
E' . The physical energies are obtained from E n = E' — E vac . The mass 



spectrum is then given by M n = y E 2 — P 2 . The low lying states of the 
massive Schwinger model are a vector state and next a scalar state |T^, KM. 
The behaviour in the ultra- relativistic region m/g — > has been computed 
in chiral perturbation theory up to second order by Adam |Tj| and Vary et 
al. [17|. The vector particle mass behaves as 



A = 2e 7 « 3.5621, B « 5.4807, C » 2.0933, (2) 

where M^ ^ = gj\pH is the boson mass of the massless Schwinger model. 
The mass of the vector boson in the chiral region is shown in Fig. [1]. The 
vector particle is almost entirely a fermion-antifermion (qq) bound state. In 
our lattice results presented in Fig. [1] only this qq sector has been taken into 



account. We find good agreement with chiral perturbation theory [|Tj| [L7j and 
with finite lattice results by Hamer et al.[]14||. Fitting our data to a second 
order polynomial in m/g yields the coefficients A = 0.56, B = 1.74 and C = 
0.22. Results for the binding energy M v — 2m of the vector boson in the non- 
relativistic region are shown in Fig. [2]. We compare our results with finite 



lattice results [0 and with results by Mo and Perry |y| using light-cone 
quantization with an expansion in orthogonal functions. The corresponding 
results for the binding energy M s —2m of the scalar boson are shown in Fig. [3]. 
Our results correspond to the qq sector. For small masses, i.e. m/g w 1/2 
to 1/8, the qqqq sector becomes important. Thus our results have a larger 
error in that region. Results on the scalar mass in the ultra-relativistic region 
m/g << 1 taking into account the qqqq sector are shown in Figs. [4,5]. 
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3.2 Dependence on #-angle 

The massive Schwinger model has #-vacua and one can study its ^-action. 
The ^-physics of the massive Schwinger model on a circle has been discussed 
by Manton ||. The wavefunction is invariant under local gauge transforma- 
tions. The #-angle is a parameter, which characterizes the behavior of the 
wavefunction under global gauge transformations *f?[A] — ► e ^[A], where 
n = 0, ±1, ±2, • • •. One can also transform Hilbert space such that this 
common phase factor disappears and the wavefunction becomes completely 
invariant under local as well as global gauge transformations. This corre- 
sponds to introducing a ^-action S[A, 9] = S[A] — 9q[A], where q[A] is the 
topological charge, or an equivalent ^-Hamiltonian. The mass spectrum of 
low-lying states as a function of the 0-angle is shown in Figs. [4,5]. Fig. [4] 
shows results (raw data), where the usual vacuum energy E vac {9 = 0) has 
been subtracted in all ^-sectors, which is unphysical, e.g., for 9 = it. It shows, 
however, the behavior of the vector mass for 9 = tt when increasing the lat- 
tice from N = 6 to 24. One observes a cusp in the curves, becoming more 
pronounced when N increases. This behavior is due to level crossings with 
higher lying states (not shown in the Fig.). Fig. [5] shows the same results 
when the vacuum energy E vac {9) has been subtracted in the corresponding 
^-sector. One observes agreement with first order chiral perturbation theory 
in the ultra-relativistic regime (m/g < 0.04 for N = 6 and m/g < 0.01 for 
iV = 24). The slope of those curves is related to the expectation value of 
the chiral condensate. The Feynman-Hellmann theorem relates the fermion 
condensate < ipip > to the derivative of the vacuum energy with respect to 
the fermion mass. Adam JL5|| has computed the condensate in lowest order 



of chiral perturbation theory, 



a = — ^ 0.2835, b w 0.7825, c w 0.7163. (3) 
2tt 



Eqs.(gD imply 



t , d „ r 

2n am 



(4) 

m=0 

which in analogy to the Feynman-Hellmann theorem relates the condensate 
to the vector mass. When we extract the slope dM v /dm from Fig. [5], we 
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obtain < ipip > /g — 0.16 cos(#), compared to the exact result of the massless 
Schwinger model given by 

( : cos(fl) ^ 0.1599 cos(0). (5) 



The ^-dependence and global gauge transformations have been treated in 
standard lattice gauge theory based on compact link variables U^(x) by 



Hamer et al. [18|. Our results are in agreement with theirs. 



3.3 Distribution functions 

In the Hamiltonian approach it is easy to compute the wave function of a low- 
lying state. From the wave function one can obtain information on its parton 
structure, i.e., the number of partons and their momentum distribution. The 
distribution function of the vector boson is given by 

f(x B ) =< | l - [b\b k + 44] I >, (6) 

where xb = k/P is the fraction of momentum of the vector boson carried by 
the parton, i.e., fermion. In Fig. [6] we display our lattice results for m/g = 
to 0.28. In the massless limit the distribution function has the shape of a 
box. Our results are compared with variational calculations using the infinite- 



momentum frame by Bergknoff [IIJ and also with those by Mo and Perry 
[ |l3|j using the light-cone. The most sensitive region is the ultra-relativistic 
region. We find agreement in shape with Mo and Perry's results and very 
good agreement with Bergknoff's results as shown in Fig. [7]. The agreement 
with Mo and Perry's results becomes even better than the agreement with 
Bergknoff if their distribution function is normalized to one. In Fig. [8] we 
also display the convergence behavior of the vector distribution function for 
m/g = 32, i.e., in the non-relativistic region. It can be seen that a reasonably 
large lattice (N ~ 200) is needed to resolve this function (which has a sharp 
peak at xb = 1/2 when m/g — > oo). 



4 Summary 

In this work we have applied a Hamiltonian lattice approach in a fast moving 
reference frame to study the massive Schwinger model. We find that the 
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method works well for the computation of the low-lying mass spectrum, also 
in the presence of the ^-action, and distribution functions. Here we have 
investigated only 9 = 0, n. It is interesting to further study the dependence 
on other ^-angles and higher-lying states. More ^-angles can be studied, e.g. 
by adding one negative momentum state to the basis. Those investigations 
are feasible, but would require a higher computational effort. 
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Figure Caption 

Fig.l The dimensionless mass of the vector boson versus m/g in the chiral 
region. Comparison with Hamer et al. |T3|] and chiral perturbation 



theory 



Fig. 2 The dimensionless binding energy of the vector boson versus 
log 2 (m/g) in the non-relativistic region. Comparison with Hamer et 
al. and Mo and Perry @. 

Fig. 3 Same as Fig. [2], but for scalar boson. 

Fig. 4 Dependence of the dimensionless mass of vector and scalar particle on 
the 6*-angle before proper subtraction of vacuum energy (raw data). 



Comparison with chiral perturbation theory JTq, 17 



Fig. 5 Dependence of the dimensionless mass of vector and scalar particle 
on the #-angle after subtraction of ^-vacuum energy. Comparison with 



chiral perturbation theory |L6|, [17 



Fig. 6 Distribution function of the vector boson for different m/g in the chiral 
region. 

Fig. 7 Distribution function of the vector boson for m/g = 0.1/ ^/tt. Com- 
parison with results by Mo and Perry Jl3| and Bergknoff |10|| . 

Fig. 8 Distribution function of the vector boson for m/g = 32 in the non- 
relativistic region. 



8 



The Mass of the Vector Boson 

In the Chiral Limit 
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The Low-Lying Mass Spectrum- Raw Data 
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Distribution Function of the Vector Boson 

Comparison with Variational Computation 







™ = °" 1 =0.0564.. 

9 ^ 




— 


/ ^r>^ 

/ xO 
is 

I 




V VA A \ 

<\\ 
1 


< 


1 


Bergknoff 

Mo and Perry 

O This work 






1 i 1 


I 1 1 


1 1 1 



0.0 0.2 0.4 0.6 0.8 1.0 

X B 



